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In this paper we have considered the Universe to be filled with Modified Gas and the Cosmo- 
logical Constant A to be time-dependent with or without the Gravitational Constant G to be 

. 2 

time-dependent. We have considered various phenomenological models for A , viz., A « p,A oc ^ 
f~ ^ ■ and A oc |. Using these models it is possible to show the accelerated expansion of the Universe at 

the present epoch. Also we have shown the natures of G and A over the total age of the Universe. 
, Using the statefinder parameters we have shown the diagramamtical representation of the evolution 

of the Universe starting from radiation era to ACDM model. 
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I. INTRODUCTION 

There are two parameters, the cosmological constant A and the gravitational constant G, present in 
Einstein's field equations. The Newtonian constant of gravitation G plays the role if a coupling constant 



between geometry and matter in the Einstein's field equations. In an evolving Universe, it appears natural 



to look at this "constant" as a function of time. Numerous suggestions based on different arguments 
have been proposed in the past few decades in which G varies with time [1]. Dirac [2] proposed a theory 
with variable G motivated by the occurrence of large numbers discovered by Weyl, Eddington and Dirac himself. 

It is widely believed that the value of A was large during the early stages of evolution and strongly influenced 
its expansion, whereas its present value is incredibly small [3] . Several authors [4] have advocated a variable A 
in the framework of Einstein's theory to account for this fact. A as a function of time has also been considered 
■ in various variable G theories in different contexts [5]. For these variations, the energy- momentum tensor of 
1^ \ matter leaves the form of the Einstein's field equations unchanged. 

o 

In attempt to modify the General Theory of Relativity, Al-Rawaf and Taha [6] related the cosmological 
constant to the Ricci Scalar 1Z. This is written as a built-in-cosmological constant, i.e., A oc 1Z. Since the 
Ricci Scalar contains a term of the form -, one adopts this variation for A. We parameterized this as A oc -. 
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Similarly, we have chosen another two forms for A : A oc p and A oc % ; where p is the energy density. 

X : 

Recent observations of the luminosity of type la Supernovae indicate [7, 8] an accelerated expansion of the 
Universe and lead to the search for a new type of matter which violates the strong energy condition i.e., p+3p < 0. 
The matter consent responsible for such a condition to be satisfied at a certain stage of evolution of the universe 
is referred to as a dark energy. There are different candidates to play the role of the dark energy. The type of 
dark energy represented by a scalar field is often called Quintessence. The simplest candidate for dark energy 
is Cosmological Constant A. In particular one can try another type of dark energy the so-called Chaplygin gas 
which obeys an equation of state like [9] p = —B/p, (B > 0), where p and p are respectively the pressure and 
energy density. Subsequently the above equation was modified to the form [10] p = —B/p a ,0 < a < 1. There 
are some works on Modified Chapligin Gas obeying equation of state [11, 12] 

p = Ap-B/p a ,(A>0) (1) 
In this work we have considered a cosmological model for the cosmological constant of the forms : A oc -, 
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A oc % and A oc p in presence of Modified Chaplygin Gas, with or without the variation of Gravitational 
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Constant G. 

In 2003, V. Sahni et al [13] introduced a pair of parameters {r, s}, called statefinder parameters. The 
trajectories in the {r, s} plane corresponding to different cosmological models demonstrate qualitatively different 
behaviour. The above statefinder diagnostic pair has the following form: 

r - 1 

and s = — -, p- (2) 



aH* 3(,-I) 

where H (= ^) and q (= — |§) are the Hubble parameter and the deceleration parameter respectively. 
The new feature of the statefinder is that it involves the third derivative of the cosmological radius. These 
parameters are dimensionless and allow us to characterize the properties of dark energy. Trajectories in the 
{r, s} plane corresponding to different cosmological models, for example ACDM model diagrams correspond to 
the fixed point s — 0, r = 1. 

The paper is organized as follows: Section II deals with Einstein field equations in presence of dynamic 
cosmological constant A. In sections III and IV we have considered different A-dependent models without and 
with variable gravitational constant G respectively and shown various stages of the evolution of the Universe 
for these models using statefinder parameters. We have taken some particular values of the proportionality 
constants for the graphical representations. We have discussed this particularization of the constants and their 
physical consequences in section V. 



II. EINSTEIN FIELD EQUATIONS AND DYNAMIC COSMOLOGICAL CONSTANT 

We consider the spherically symmetric FRW metric, 

dr 2 



ds A =dt A -a 2 (t) 



kr 2 



+ r 2 (d9 2 + sin 2 6d(f> 2 ) 



(3) 



where a(t) is the scale factor and k is the curvature scalar with values 0, 1 and —1 for respectively flat, closed 
and open models of the Universe. The Einstein field equations for a spatially flat Universe (i.e., taking k = 0) 
with a time-dependent cosmological constant A(t) are given by (choosing c = 1), 



and 



3^J = 8irGp + A(t) (4) 
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2- + ^ = -87rG P + A(t) (5) 



where p and p are the energy density and isotropic presure respectively. 

Let us choose modified chaplygin gas with equation of state given by equation (1). 

Here, we consider the phenomenological models for A(t) of the forms A oc p, A oc % and A oc -. 
First we will consider G to be constant and try to find out the solutions for density p and the scale 
factor a(t) and hence study the cosmological models in terms of the statefinder parameters r, s. Secondly 
we will consider G to be variable as well and study the various phases of the Universe represented by the models. 



III. MODELS KEEPING G CONSTANT AND A VARIABLE 

Taking G to be constant and A to be time dependent, the energy conservation equation is, 

P + i a {p + P) = -^G 
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Fig. 1 shows the variation of s against r for different values of /3i = 1, § respectively and fora = 1, A = 1/3 and 

8 7T G = 1. 



A. Model with A oc p 

Here we consider 

A = /?i p 

where (3\ is a constant. 

Equation (7) together with equations (1) and (6) yield the solution for p to be, 



(7) 



B 



C 



I I 1 I A 24irG(l + A)(l + c) I 

+ a 8*<2+Pi / 

where G is an arbitrary constant. 

Substituting equation (7) and (8) in equation (4), we get the solution for the scale factor a(t) as, 



a hh y/8irG + ]h 2 F 1 [/ 2 , / 2 , 1 + f 2 , 



,/iB 



G(l + A) 



] = 4V3(1 + A)GirC f2 t 



(8) 



(9) 



24ttG(1 + ,4)(1 + o 



where /i = an d /2 = 2(i+a) • Hence, for small values of a(t), we have, p ~ 

which is very large and the equation of state (1) reduces to p ~ Ap. Again for large values of a(t), we get 
i i 

p ~ (iTa) 1+ ° an< ^ V — — (iTa) 1+ ° ' ^' e '' P — ~P which coincides with the result obtained for MCG with 
/?i=0 [12]. 

Using equations (2), (4) and (6) we get, 

36ttG(1 + y)[8TTG{A(l + a) - ya} - fa] 



1 + 



(87TG + /3!) 2 

2 furthei 

Therefore for acceleration q < 



8ttG(1 + j/)[87rG{A(l + a) - ya} - ft] 
* = (8.G + / 3l )(8 7 rG,-/3 1 ) (10) 



where y = jjj which can be further reduced to a single relation between r and s. Now q — — — s^G(i+3y) 



y < j§^q — |- Also for the present epoch q — 



«7P ~~ 2(87rG+,3i) 
2 ^ f = 3^8^G 



1/ 3/3i _ 9 n 
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Fig.2 

Fig. 2 shows the variation of A against a(t) for different values of f}\ = 1, § respectively and for a = 1, A = 1/3, 8 n G 
= 1, B = 1, C = 1 . 



If we assume that the present Universe is dust filled, we have y = 0, i.e., /?i = 16 g G - Taking 87rG = 1 
we get the best fit value to be f3\ = |, which gives r = 1 ( choosing A = i, a = 1 ) for the present 
Universe. That means the dark energy responsible for the the present acceleration is nothing but A. Also 
A = 1, a = 1 and /?i = | give r = 2.16 for the present time. For this case /?i > 3 gives non-feasible 
solutions in the sense that the present values of y, i.e., £ becomes too large. For f3\ = 1, we get the 

present value of y to be |, but again r = 1. In cither of the above cases we get accelerating expansion 
of the Universe. These can be represented diagrammatically in the r, s plane. This is shown in figure 1 
(taking A = |,a = l,/?i = l,|,87rG = 1 and A = l,a = l,/3i = 1, |,87rG = 1). Figure 1 represents the 
evolution of the Universe starting from radiaton era to ACDM model. Here we get a discontinuity at (3\ = —8irG. 

Again for this model 

( B C \ ^ 

A = M TT~T + Eimssis ( n ) 

Variation of A(t) against a(t) is shown in figure 2 for different choices of which represents that regardless 
the values of fa, A(i), i.e., the effect of the cosmological constant decreases with time. 



B. Model with A oc H 2 

Choosing 

A(t) = /? 2 # 2 (12) 
where /? 2 is a constant and proceeding as above, we obtain the solutions for p, a(t), A as, 

(B C \ 

TTA + a (3-/32)(l+A)(l+a) J ^ 13 ) 

a hh 2 F 1 [f 2 , f 2 , 1 + h, - C V +A) \ = V^GV^lh (1 + A)C^ t (14) 



where h = (3 - /3 2 )(1 + + a) and / 2 = 



8irG(3 2 ( B C , 
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Fig. 3 shows the variation of s against r for different values of /? 2 = 1, 2 and for a = 1, A = 1/3, 8 it G = 1. Fig. 4 
shows the variation of A against a(t) for different values of /?2 = 1, 2 and for a = 1, A = 1/3, 8 n G = 1. 



Here for fa < 3 we can check the consistency of the result by showing p ~ at small values of a(t) and 
p = — p for large values of a(t). But if we take fa > 3 we get opposite results which contradict our previous 
notions of the nature of the equation of state (1). Again for fa = 3, we get only ACDM point ,i.e., we get a 
discontinuity. Therefore, we restrict our choice for fa in this case to be fa < 3. 

Computing the state-finder parameters given by equation (2), we get the equations for r and s to be, 

(3- fa){l+y)[{A{l + a)-ya}(3-fa) -fa\ _ (3 - fa)(l + y)][{A(l + a) - ya}(3 - fa) - fa] 

+ 2 H(3-fa)y-fa} 

(16) 

(where y = ^), which can still be resolved into a single relation and can be plotted in the r,s plane. 
Here q = |[(3 — fa)y — {fa — 1)]. Hence the Universe will accelerate if q < y < fz^j- Again 

for the present Universe g = — ^ ^> y = f 2 ^- Assuming the present Universe to be dust dominated, 

i.e., y — we get the best fit value for fa to be 2. Taking A = = l,8nG = I, fa = 2 and y = 0( 

1. e. dust dominated present Universe ) we get the present value to be r = 1/3, also the same values with 
fa = 1 gives the present values to be y = — i, r = | .This is shown in figure 3 (A = h,a = 1, fa = 1 and 

2, 87rG = 1) , which explains the evolution of the Universe from radiation era to ACDM model. Again variation of 
A against time is shown in figure 4, where we can see that A decreases with time for whatever the value of fa be. 



C. Model with A oc ^ 

a 

Taking 



A = fa- (17) 
a 

(where fa is a constant), and proceeding as above we get a relation for p as, 

Unlike the previous two cases here we get a far more restricted solution. Here the only choice of fa for which 
we get the feasible solution satisfying p ~ Ap for small values of a(t) and p~ —p for large values of a(t) is 

fe< (l + A)(l + ^) ° r P3>3 (19) 
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Fig. 5 shows the variation of s against r for different values of A = i, a = 1, SttG = 1,03 = |, 4 and | 



Again since g = — -^jp = — 2 = 4 ( 3 _^j"_H^ = 4 ( 3 -fe^H^ ; /? 3 > 3 implies q < without even violating the 

energy-condition p + 3p > 0. Although /3 3 < ( 1+ ^( 1+a ) causes the acceleration of the Universe violating the 

energy-condition. Taking q = — | for the present epoch, we obtain y = 7§z^y- Hence the present epoch is dust 

filled if (33 — 4 and thus giving the present value of r to be — ^ for A = 3, a = 1, 87rG = l.On using relation 
(19), p and therefore a, A cannot be expressed in an open form. We can rather derive a solution for A in terms 
of p, p as, 

A=|^(P + 3 P ) (20) 

P3 - o 



Using equations (2)we get the statefinder parameters as, 



_ i (l + y)(/3 3 -3)[/3 3 + (/?3 + 6) a; ] 2(l + j / )(/3 3 -3)[/3 3 + (/3 3 + 6)x] 

(/?3 + /3 3 x- 2)(/3 3 +/3 3 2/-2) [ / 33 + ( /3 3 + 6 )y][ /3 3 + / 3 3a; _2] 1 j 



where y = ^ and a; = i.e., x = A(l + a) — ya [from equation 



(i)]. 



Eliminating y between the equations (21), we get a single relation of r and s, which can be represented 
diagrammatically in the r,s plane (figure 5). Here we have taken A = = l,8nG = l,/3 3 = ^,4 and |, 
combining two cases. Taking /3 3 = |, | we can explain the evolution of the Universe starting from - = — | to 

ACDM model and /3 3 = 4 explains the evolution of the Universe starting from radiation era to y = — | , as seen 
from the expression for q. Considering the present epoch to be dust-dominated, the present value of r is given 
for /3 3 = 4 to be — |. As follows, the former two cases cannot give the present value of r, as y = > — | for the 
present epoch. Here we have an infinite discontinuity at ^ = — |, i.e., when p + 3p = 0. Also since we do not 
get a closed from of p here, it is difficult to plot A against the scale factor a(t) . 



IV. MODELS WITH A AND G BOTH VARIABLE 

Now we consider G as well as A to be variable. With this the equations (3), (4), (5) yield the conservation 
laws as, 

p + 3H(p + p) = (22) 

and 

A + 8nGp = (23) 
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Fig. 6 shows the variation of s against r for different values of 71 = 1,3 and 3.5 and A = |, a = 1, B = 1, Ci = 1. 



Now we study the various phases of the Universe represented by these models. 
Equation (22) together with equation (1) yield the solution for p as, 



1 „3fl+v41fl+«'> 

(24) 



1 + A a 3(l+^)(l+a) 

where C is an arbitrary constant. This result is consistent with the results already obtained [12]. 



Here we consider 

where 71 is a constant. 

Equation (22), (23) and (25) give, 



A. Model with A oc p 



A = 71 P (25) 



8tt 

where C\ is a constant and p is given by equation (24). 
Using equations (2), (4), (23) and (26), we get, 

G = Ci + 2i ^±«) log( _B_ ) 

r _ 1 j_ 9(l+y)[87rG{A(l+ a )-y a }- 7l (l+ y )] 



G = Cl -^\ogp (26) 



2(8^G+ 7 l) ~ " (27) 

_ (l + y)l87TG{A(l+ a )-ya}- 7l (l+y)] 
* _ (877Gy- 7 l) 

where y = -. 

Equation (27) cannot be resolved to get a single relation between r and s, rather we obtain a parametric 
relation between the same with y — - as the parameter. This can be represented diagrammatically in the 
r,s plane, which is shown in figure 6 taking ji = 1,3 and 3.5 and A = |,a = 1,B = l,C\ = 1. Now 
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Fig. 7 shows the variation of A against a(t) for different values of 71 = 1, 3, 3.5 and for a = 1, A = 1/3, C\ = 1. Fig. 8 
shows the variation of G against a(t) for different values of 71 = 1, 3, 3.5 and for a — 1, A = 1/3, 

Ci = 1,B = 1, C = 1. 



Q = 47TG s^G+'y 1 71 • Taking into account that q = — | for the present epoch, wc get y = — |. Therefore, 
for the present dust-dominated era y — and 71 = 16 g G . hence for the This models represents the Universe 
starting from the radiation era to ACDM model. Again figure 7 represents the variation of A against the scale 
factor a(t) with 71 = 1, 3, 3.5 and figure 8 represents the variation G against the scale factor a(t). These figures 
show that for this particular phenomcnological model of A, G starting from very low initial value increases 
largely and becomes constant after a certain period of time, whereas A starting from a very large decreases 
largely to reach a very low value and becomes constant. 



B. Model with A oc H 2 

We consider 

A = l2 H 2 (28) 

Proceeding as above we get, 

A = 8ttG—^—p (29) 

3-72 

where 72 is a constant. 

Solving equation (22), (23) and (29) we get, 

G=-^ (30) 
P 3 

where C2 is a constant. 

Using equations (2), (4) and (23), we find the state-finder parameters as, 

r = 1 H (1 + I/)(3 - 72)[3{A(1 + a)- ya} - (1 + 37)72] (31) 

(1 + y)(3 - 72 )[3{A(1 + a) - ya} - (1 + y) l2 ] 
(3 - 72)2/ - 72 
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Fig. 9 shows the variation of s against r for different values of 72 = 1 and 3.5 and j4=i,a = l,S = l,(72 = l, C = 1. 
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Fig. 11 



Fig. 10 shows the variation of G against a(t) for different values of 71 = 1, 3.5 and for a — 1, A = 1/3, 
Ci = 1, C = 1, B — 1. Fig. 8 shows the variation of A against a(t) for different values of 71 = 1, 3,3.5 and for a 

1, A = 1/3, C 2 = 1,5 = 1, C = 1. 



where y = -. 

Now g = i[(3 — 72)2/ — (72 — 1)]]. These equations can further be resolved into a single relation of r and 
s, which can be plotted diagrammatically in the r, s plane. Here we get a discontinuity at 72 = 3. We have 
plotted these values in the r, s plane taking 72 = 1 and 3.5 in figure 9 (A = |, a = 1). This case explains the 
present acceleration of the Universe, starting from radiation era to ACDM model. 

Also figures 10 and 11 show respectively the variation of G and A against the scale factor for the same values 
of the constants. Here also like the previous case G starting from a very low initial value increases largely and 
then continues to be constant near unity. On the other hand A starting from a large value decreases largely 
and continues to be constant after a certain period of time. 



C. Model with A oc - 

a 

Here we consider 



A = 73- (33) 
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Fig. 12 shows the variation of s against r for different values of 72 = 2 and 3.5 and 
A= |,a = 1,B = 1,C 3 = 1, C = 1. 



where 73 is a constant. 



Using equation (33) in equations (4) and (5), we get, 

3-73 

Also, G can be solved to be, 

G = C 3 [p 2 -^(i+A) {2 - 73(1 + A) + ^a+^) + (i+ Q )(2- T3 (i+A))>]- (35) 

Using equations (2), (4), (22), (23), we find, 

_ . , (1 + y)(3 - 73)[6{A(1 + q) - ya} + 73 (1 + y)] 

+ [2 -73(1 + W 



2(1 + y)(3 - 73)[6{A(1 + a) - ya} + 73 (1 + y)} 
3[2-7 3 (l + y)][73 + (73 + 6)2/ 

where y — 2 and C 3 is a constant. Equations (36) and (37) can further be resolved to get one single relation 
between r and s and plotted diagrammatically taking 73 = 2 and 3.5 (figure 12). Since deceleration parameter 
q = —-^jj>z = — 73 X H i = ( 3 ** 2) H-i 1 is negative in the present epoch, we get 3 — 73 < 0, i.e., 73 > 3. Also for 
73 = 3 we get discontinuity. Both the models represent the phases of the Universe starting from radiation era 
to ACDM model. Again G and A can be plotted against a (figures 13 and 14 respectively) . Unlike the previous 
cases this model an opposite nature of G and A, as G decreases with time and A increases with time. 



V. DISCUSSION 



Here we have considered three phenomenological models of A, with or without keeping G to be constant. 
Keeping G constant we always get accelerated expansion of the Universe. For the first case, i.e., A oc p 
or more precisely, A = /3\p, for particular choices of the constants we get that the dark energy responsible 
for the present acceleration is nothing but A. Also the density parameter of the Universe for this case 
is given by, flf^l — 8 ^-f = 87r 8 J^ /3l and the vacuum density parameter is Sl^ 1 = ^jj? = S7T § 1 +(3l , so that 
Qtotai = + ^A = + = 1- Also for A oc H 2 , i.e., A = /?2 H 2 , the density parameter and vacuum 
density parameter are given by, il!^ = ^p 2 - and fi^ 1 = ^ respectively, so that Qtotal = + ^a = 1- Again 
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Fig. 13 shows the variation of G against a(t) for different values of 71 = 2, 3.5 and for a — 1, A = 1/3, 
C3 = 1, C = 1, B — 1. Fig. 14 shows the variation of A against a(t) for different values of 71 = 2, 3, 3.5 and for a 

1, A = 1/3, C 3 = 1,B = 1, C = 1. 



for A cx ^ or A = /3 3 ^ , we have the corresponding parameters as, 0^ = g^^^^E) 1 ^a* = 3(2- 3 fe-fe£) an< ^ 
^totai = 1- Now ^tota/ = fi{^ + fi^ 3 = 1 for all the models. Also we can compare these models by taking, 
n^l = fi^, so that /?2 = rsnG+M ' ^ ow we wom d like to take into account the present values of the density 
parameter and vacuum parameter obtained by the recent measurements. Considering f2 m o = 0.33 ± .035, we 
calculate the present values of the proportional constants to be 1.7397if < 0i < 2.3898X, 1.905 < /3 2 ° < 2.115 
and 3.7937 < /3 3 ° < 4.2099, where K — 8ttG and G is the present value of the gravitational constant. Thus 
we get the value of P3 to be lesser than the previous works [5, 14]. Again considering G to be time-dependent, 
we get the same values of the parameters as that with G constant, i.e., the ranges of 7i°,72°7 73° are same as 
that of (3\ , (32°, /?3° respectively. Here also we get cosmic acceleration and the nature of variation G and A 
as well. We get two different cases regarding the variation of G and A. For the first two cases we see that G 
increases and A decreases with time, whereas for the third case G decreases and A increases with time. In all 
the cases the values become constant after a certain period of time, i.e., the present day values of G and A are 
constants. Thus these models with the phenomenological laws give us some interesting features of the cosmic 
acceleration and some modified values of the parameters. Also we get the natures of the Cosmological Constant 
and the Gravitational Constant over the total age of the Universe. We can also make use of the statefindcr 
parameters to show the evolution of the Universe starting from radiation era to ACDM model. 
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